Preprint typeset in JHEP style. - HYPER VERSION 

Deviations From Newton's Law in 
Supersymmetric Large Extra Dimensions 



p. Callin^ and CP. Burgess^'^ 

^ Institute of Physics, University of Oslo, N-0316 Oslo, Norway. 
^ Department of Physics and Astronomy, McMaster University, 

1280 Main Street West, Hamilton, Ontario, Canada, L8S ^Ml. 
^ Perimeter Institute, 31 Caroline Street North, 

Waterloo, Ontario, Canada. 



E-mail: a.p.callin@fys.uio.no, cburgess@perimeterinstitute.ca 



Abstract: Deviations from Newton's Inverse-Squared Law at the micron length 
scale are smoking-gun signals for models containing Supersymmetric Large Extra 
Dimensions (SLEDs), which have been proposed as approaches for resolving the 
Cosmological Constant Problem. Just like their non-supersymmetric counterparts, 
SLED models predict gravity to deviate from the inverse- square law because of the 
advent of new dimensions at sub-millimeter scales. However SLED models differ 
from their non-supersymmetric counterparts in three important ways: {%) the size 
of the extra dimensions is fixed by the observed value of the Dark Energy density, 
making it impossible to shorten the range over which new deviations from New- 
ton's law must be seen; {ii) supersymmetry predicts there to be more fields in the 
extra dimensions than just gravity, implying different types of couplings to matter 
and the possibility of repulsive as well as attractive interactions; and {Hi) the same 
mechanism which is purported to keep the cosmological constant naturally small also 
keeps the extra-dimensional moduli effectively massless, leading to deviations from 
General Relativity in the far infrared of the scalar-tensor form. We here explore the 
deviations from Newton's Law which are predicted over micron distances, and show 
the ways in which they differ and resemble those in the non-supersymmetric case. 
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1. Introduction 

The recent discovery of a non-vanishing Dark Energy density [|T| throws into sharp 
relief the general theoretical befuddlement about how to predict reliably the gravi- 
tational response to the energy of the vacuum p|. In a nutshell, the difficulty arises 
because in four dimensions gravity responds to the energy density of the vacuum, 
p, as if it were a cosmological constant. But p generically receives its largest contri- 
butions from the quantum zero-point energy of the shortest-wavelength modes, for 
instance with modes having wavelength £ contributing in four dimensions an amount 
of order 5p ~ As such, the theoretical prediction is sensitive to the most micro- 
scopic of a theory's details, and since any theory describing the known particle types 
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has I ^ larger than 100 GeV, the theoretical prediction is generically many orders of 
magnitude larger than the observed value: ^ 

The lie 

Pobs = with 3 X 10"^ eV. (1.1) 

a 

1.1 Naturalness Issues 

Until the observational result was found to be nonzero, the theoretical goal had been 
to identify a theory which is both experimentally successful and technically natural 
in that sense that the contributions to the vacuum energy for some reason cancel 
once summed over the relevant particle content. Now that we know the result is 
nonzero the bar has been raised, and a technically natural solution would instead 
require a cancellation only for those modes having scales i < a. Although this kind 
of naturalness is not absolutely required by fundamental principles, it amounts to the 
requirement that the small hierarchy can be understood within an effective theory 
of microscopic physics, regardless of the scales for which this theory is formulated. 
Technical naturalness is a conservative requirement because we know that it applies 
for all of the many other hierarchies of scale for which we have solid evidence, in all 
branches of physics. 

Yet it is hard to modify physics at scales above ^ ~ a to get a small enough 
vacuum energy without also running into conflict with the many experimental mea- 
surements which are available for such scales. More recently, despair at making 
progress along these lines has led many to abandon the criterion of technical natural- 
ness altogether and to turn to anthropic ideas along the lines proposed by Weinberg 
more than ten years ago 0, This point of view has gained additional momentum 
from the observation that anthropic reasoning might be required to make predictions 
within string theory, given the enormous number of vacua which arise there . This 
line of argument has emboldened others to abandon using technical naturalness as a 
criterion for understanding the electro- weak hierarchy 0. 

On the other hand, potential progress has also been made in obtaining a small 
gravitational response to the vacuum energy without abandoning technical natu- 
ralness, within the framework of Supersymmetric Large Extra Dimensions (SLED) 
[0, ^ . The idea behind this proposal is that it is indeed possible to modify physics 
at scales ^ < a in such a way as only to modify gravitational interactions without 
changing non-gravitational physics. As a result, although the vacuum energy may 
^From here on we adopt units for which h = c — 1. 



2 



be large it may be possible to modify how this vacuum energy gravitates for scales 
i < a without also ruining the agreement between theory and the vast number of 
non-gravitational experiments. The framework within which this can be done is that 



of large extra dimensions |T0[, according to which two internal dimensions are taken 
to be as large as they can possibly be without coming into conflict with observations. 
Interestingly enough, the largest radius which is possible for such extra dimensions 



turns out to be r ~ a, as imposed by tests of Newton's Law at short distances |]TT], |T^]. 

It is proposed that within such models supersymmetry in the extra dimensions 
can play the crucial role of cancelling the contributions to the vacuum energy from 
all modes having wavelengths £ < r ~ a. There is currently considerable activity de- 
voted to understanding what the gravitational response to the vacuum energy would 
look like within this framework, in order to see whether this kind of modification 
really provides a technically natural description of the dark energy J^, ||, H, |13|, |14 



1.2 The Smoking Gun 

A particularly attractive feature of the SLED framework is that it is very predictive, 
and so eminently falsifiable. This predictiveness ultimately comes because technical 
naturalness requires the proposal to modify physics at length scales £ < a (or energies 
above a^^ ~ 10^'^ eV) — a range of scales to which we have a great deal of experi- 
mental access. In particular SLED models add to the observational implications for 
the Large Hadron Collider (LHC) of the earlier large-extra-dimensional models [|1^ 



as well as predicting new features to do with supersymmetry |T7|, |T8[. As such, 
SLED models provide a unique framework which links the technical naturalness of 
the vacuum energy with observational implications at high-energy accelerators. 

Given that its raison d'etre is modifying the gravitational response of the vac- 
uum, it is the modifications to gravity which provide the most definitive signatures 
for the SLED proposal. These come in two types, corresponding to deviations from 
General Relativity over both sub-millimeter and cosmological distance scales. Over 
very large distances the model behaves like a scalar-tensor theory, with the light 
scalar (or scalars) describing the dynamics of the moduli of the extra dimensions. 
Unlike for non-supersymmetric models having large extra dimensions, within SLED 
models there is an understanding of why these modes remain essentially massless 
in a technically natural way, which is a consequence of the same mechanism which 



is purported to protect the vacuum energy in these models [O]. This leads to im- 
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plications for the time- dependence of the Dark Energy density |T^, |20|, as well as 
model-dependent implications for long-distance tests of gravity. 



However, it is tests of Newton's Law over micron length scales which ultimately 
provide the smoking gun for SLED models. This is because SLED models rely for 
their success on the existence of two large extra dimensions, and on the coincidence 
between the size, r, of the extra dimensions and the observed scale, a, in the observed 
Dark Energy density. This means that Newton's law must be violated at length 
scales. A, which are of order the Kaluza Klein (KK) masses for the large extra 



dimensions: A ~ r/27r ~ a/27r, which works out to be around a micron in size [|T7 
Since the size of these dimensions cannot be shrunk without making the observed 
vacuum energy too large, there is no way to escape the implication that Newton's 
law must change once tested over these distances. 

Since these deviations are predicted over distances which are not too far from 
experimental reach, we can hope that the theory can be definitively tested within the 
not-to-distant future. For these purposes it is important to have precise predictions 
for what kinds of deviations should be expected. It is the purpose of this paper to 
provide a first calculation of the short- distance behaviour of force laws within the 
SLED picture, working within the framework of the simplest (toroidal) compacti- 
fications. In particular, we compute the forces which should be expected between 
bodies which are separated in the visible dimensions but are not displaced in the ex- 
tra dimensions, due to the mediation of the various bosons which propagate within 
the extra dimensions. As might be expected, there are more of these bosons than 
arise in non-supersymmetric models due to the additional particle content which 
supersymmetry in the bulk requires. Some of these additional bosons can mediate 
spin-dependent forces, and the detection of such could provide a way to distinguish 
between the supersymmetric from non-supersymmetric options. 

We organize our results as follows. The next section, §2, briefly describes the 
bosonic field content of 6D supergravity theories, and then describes the quantization 
of their linearized fluctuations about flat space. The exchange of these fluctuations is 
then used in §3 to compute the interaction between slowly- moving, localized classical 
sources, and identify their interaction potential energy. We also identify what kinds 
of charges such sources can carry in the static limit, and follow how the interaction 
energy depends on these charges. Finally, we close in §4 with a concluding discussion. 
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2. Forces Mediated by 6D Bosons 

In this section we describe the bosonic fields which make up the bulk sector in 
models with supersymmetric large extra dimensions. We then identify how these 
bulk particles can couple to brane degrees of freedom, and compute the forces which 
are obtained when they are exchanged by two sources which are both localized on 
the same brane. 



2.1 6D Supergravity 

The field content of practical interest for SLED models is the bosonic sector of 6D 
supergravity coupled to various forms of super-matter. In (2,0) supersymmetry the 
bosonic field content of the supergravity multiplet consists of the metric, Qmni and 
a 2-form gauge potential which is subject to a six-dimensional self-duality condition. 
This is normally combined with a tensor multiplet, whose bosonic content contains 
another 2-form potential having the opposite self-duality property, and so which 
combines with the previous one into an unconstrained potential, Bmn- The other 
boson in the tensor multiplet is a scalar field, 0, known as the dilaton. The most 
common types of (2,0) matter multiplets are gauge multiplets — whose bosonic parts 
consist of the gauge potentials, A"^, of some gauge group — or hypermultiplets — 
whose bosonic part consists of scalar fields, These latter fields parameterize a 
coset space, G/ H, and so nominally describe the Goldstone modes for the symmetry 
breaking pattern G —>■ H. For some supergravities a scalar potential, f ($), can exist 
for these fields, in which case the symmetry G is explicitly broken and the fields 
become instead pseudo- Goldstone particles. Alternatively, if a subgroup of G is 
gauged, then some of these hyper-scalars can be eaten by the gauge fields through 
the usual Anderson-Higgs-Kibble mechanism. 

For a broad class of supergravities these bosons are governed by the following 
action^ 



S = d^Xy/^ 



g^'"" ( Rmn + 5a/0 9;v0 + Gam du^" dM^' 



2^2 

1 



(2.1) 



^Our metric is 'mostly plus' and we adopt Weinberg's curvature conventions which differ 
from those of Misner, Thorne and Wheeler [ p3[ only in an overall change of sign of the Riemann 
tensor. 
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where k = with M being the 6D Planck mass, which is of order 10 TeV [^]. The 
quantity Gab{^) represents a G-invariant metric on the coset space G/H, and we here 
focus our attention on the case where no scalar potential exists. Supersymmetry puts 
a number of restrictions on the form the target-space metric, Gab can take, but these 
do not play a role in what follows. Finally, F^j^ = OmA"^ — d^Al^ + c°'py A^j^^ Al^ 
is the usual gauge- covariant field strength, and the field strength, Gmnp, for the 
2- form potential, Bmn, is given by 

Gmnp = dMB^p + SnBpm + dpBMN + ck Qmnp , (2.2) 

where ^Imnp is the Chern-Simons form, given for abelian gauge groups by ^mnp = 
FmnAp + FnpAm + FpmAn- c here is a supergravity-specific constant, which can 
vanish for some of the extra-dimensional tensors appearing in the bulk supergravity. 

The equations of motion obtained from the action, ( |2.1| ), are given by Einstein's 
equation, 

1 i?^^ = dM<P dN(i> + Ga^m du^"" + \ e-2<^ GmabGn'''' - e-^ FmaFn"" 



^ 9mn 



(2.3) 



supplemented by the other field equations 



= ^ 00 + ^ e-^^ Gmnp G*^^^ + ]■ e"^ FmnF'^^ (2.4) 

= + r^,($) g^^^ 9a/$^ On^^ (2.5) 

= Vp(e-2<^G^""^) (2.6) 

= (e~<^ F^*^) + K G^'^^'Fab • (2.7) 

Here the r^^($) represent the Christoffel symbols built in the usual way from the 
target-space metric, Gabi,^)- 

As noted earlier, these expressions do not include a potential for the scalar 
fields, as would be appropriate for ungauged 6D supergravity. Even for ungauged 
supergravity it could also happen that a scalar potential is generated for fields like 
if some of the gauge field strengths, Fmn or Gmnp, are nonzero in the background. 
We do not consider here potentials coming from either sources since the success of 
the SLED mechanism relies on these potentials not being generated with a size which 
is large enough to be relevant to our present focus: the forces between macroscopic 
objects over sub-millimeter distances. 
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2.2 Toroidal Compactification 

The simplest solution to these equations is the trivial one, for which the scalars 
are constants while all gauge field-strengths vanish: — and Gmnp — 0. 

Since these assumptions ensure there is no matter stress-energy, the simplest metric 
configuration is fiat space: Qmn — Vmn- We compactify to 4 dimensions by taking 
two of these fiat directions (x"*, m = 4, 5) to fill out a compact 2-torus, 7^, of volume 
A, while the other four dimensions (x'*, fj, — 0, are large. 

For these purposes we define T2 by identifying points on the fiat plane according 

to 

{x^, x^) = (x^ + n2r2 cos 9 + nin; + n2r2 sin 6') , (2.8) 

where ni 2 are integers and 6, and r2 are the three real moduli of the 2-torus. 
Equivalently, in terms of the complex coordinate z = x"^ \ ix^ this is 

z ^ z + {n2T + ni)ri , (2.9) 

with the complex quantity r defined by r = (r2/ri) e*^ = ti + iT2, and as defined r 
takes values in the upper half-plane, T2 > 0. In terms of these quantities, the volume 
of the 2-torus becomes A — rir2 sin 9. 

To set notation, we now dimensionally reduce a free 6D scalar field on such 
a background. Consider therefore a 6D scalar field ip whose 6D field equation is^ 
(— De -l- m^) ip, propagating on a spacetime compactified to 4 dimensions on the 
above 2-torus. We assume the scalar satisfies the following boundary conditions 

i^ix^", x^ + nsrs cos 9 + niVi, x^ + n2r2 sin 9) = e^"^("i ^1+"^ ipix^", x\ x^) , (2.10) 

with < pi,2 < 1 being two real quantities. The choices pi,2 = 0, | correspond 
to taking periodic or anti-periodic boundary conditions along the torus' two cycles, 
although more general values of pi are also possible. 

The scalar field may be expanded in terms of the eigenf unctions of the 2D 
Laplacian, ^2 = 5| -|- 9|, corresponding to 

i^{x,y)^^i^p{x)up{y), (2.11) 
p 

where we denote the 4D coordinates, x'', collectively by x and the 2D coordinates, 

x^, BS y (with — x^ and — x^). The mode functions are given by plane 

^The subscript '6' here emphasizes that the d'Alembertian appearing here is the 6-diniensional 
one. 
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waves, Up{y) = A e*^'^, where the above boundary conditions imply the allowed 
extra-dimensional momenta are given by 



p ■ y = 2tt 



{ni + pi){y^ - y^cote) ^ {n2 + p2)y'^ esc 9 
ri r2 

2 , 



(2.12) 



These modes satisfy □2Mp = —p Up with 

n2 + P2- r(ni + pi 



P' 



(2.13) 



At2 

Viewed from a 4D perspective, each of the KK modes, ipp, is a 4D scalar field 
satisfying the field equation 

{-n^ + p^ + m^) ijp = , (2.14) 

and so whose 4D mass is related to the 6D mass, m, of the full 6D field ip, by Mp = 
p^ + m?. Similar decompositions also hold for free higher-spin fields compactified on 
a 2-torus, as is described in more detail (as needed) below. 

2.3 Linear Fluctuations 

The above free-field expressions are useful for the present purposes because they 
may be applied to the dynamics of small fluctuations about the flat vacuum solution 
described above. The free-field results describe the leading (linear) parts of these 
fluctuations. In particular, we wish to explore the potential energy which is set up 
by various sources through the exchange of small fluctuations in the bosonic fields 
which arise in 6D supergravity, described earlier. 

To this end we expand the action, ( |2.1| ), about the toroidal background configu- 
ration described above, working to quadratic order in the deviations of the fields 0, 

A^.j, Bmn and qmn about this background. For an appropriate choice for the 
various gauge-averaging terms (more about which below) this leads to a quadratic 
action within which fields of differing spins do not mix: 

'S'quad = + Si^ + Sa + Sb + Sh ■ (2-15) 

Each term in this action is described in more detail in the following sections. 
The Dilaton 

In the absence of background gauge field strengths, Fmn = Gmnp = 0, there is 
no potential for the dilaton and so any constant value 0o provides an equally good 
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background value. Dividing the dilaton into background and fluctuation according 
to 

(l) = (j)o + K4), (2.16) 
and expanding the action to quadratic order in the fluctuation, 0, then leads to 

S4> = -1 [ d'x v'^'^'idM^d^i) . (2.17) 



2 „ 

Since this has the same form as the example discussed above of dimensional 
reduction of a free scalar fleld on 7^ (in the special case m = 0), the results of this 
earlier section may be taken over in whole cloth. In particular we decompose in 
terms of toroidal modes as in eq. (|2.11|) , leading to the tower of KK masses given in 
eq. (|2.13| ). The momentum-space propagator for such a 6D fleld is then 

where = {fc'^, p}, with denoting the (time-like) 4-momentum in the four visible 
directions (with 4D contraction k"^ = fc'^fc^) and p representing the (Euclidean) 2- 
momentum in the compact two dimensions (with 2D contraction = p ■ p)- 

The dilaton can couple directly to brane matter, depending on the nature of the 
microscopic physics of the brane. For instance for D-branes in string theory, branes 
typically couple (in the absence of background gauge and Kalb-Ramond flelds) to 
the dilaton with strength 

St= [ df^i^/^e^^T, (2.19) 



where T is the brane tension, d is the spacetime dimension of the brane, A is a 
constant whose value depends on the kind of brane involved, and are coordinates 
on the brane world-sheet which we take to sweep out the surface x^{^) in spacetime. 
Also, 7 = det7^j, where 7^j, = gMNdf_iX^^d^x^ denotes the induced metric on the 
brane world-sheet. For some branes it can happen that A vanishes, as it does for 
instance in the case of Z?3-branes in the lOD Einstein frame. 

To linear order in the fluctuation we parameterize the dilaton interaction with 
the form 

Sint = [ d"^ ^ p , (2.20) 



where p denotes the local source density. For instance for a point particle of mass 
m on a 3-brane one might have p = X[T + m5^(^)], for some mass m. In general, 
such a coupling to the brane tension, T, can ruin the existence of the flat background 
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solution of interest here due to the back-reaction which the presence of a configuration 
of branes induces on the intervening 6D bulk fields. Indeed, an understanding of how 
the bulk adjusts is central to the issue of whether 6D supergravity provides a naturally 
small dark-energy density. In what follows we imagine that this back-reaction can 
be negligible to the linearized order to which we work, and focus on the implications 
for experiments of the dilaton couplings to localized particles. 



Goldstone Modes 

The discussion of the hyper-scalars, follows a similar route, once the fields are 
expanded into background plus fiuctuation according to 

= + K^'' . (2.21) 

The main difference is in this case the nature of the couplings to brane-bound particles 
which might be expected. 

In general, for constant background fields, c^mV^" = 0, it is always possible^ to 
rescale the fiuctuations such that Gab{f) = ^ab, in which case the quadratic action 
for <l>" becomes equivalent to N copies of the dilaton action, S"^: 

S'^ = -\j SabV^'^'dM^'^dA'. (2.22) 

This then leads to the following momentum-space propagator for <l>" 

^a.(^^P) = p^' (2-23) 

where k'^ and are defined as before. 

Because arises as a Goldstone boson its couplings are typically derivative 
couplings, which precludes there being a dilaton-like coupling of the form ( p.l9| ) 
directly to brane matter. Instead, we expect the linearized coupling to have the 
general form 

Sint = K ! d'i^ j^'du^"" , (2.24) 



where denotes the contribution of background brane fields to the conserved cur- 
rent of the symmetry for which is the Goldstone boson. 



*We assume here the target space manifold to have EucUdean signature. 
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2.4 Gauge Fields 



To quadratic order the expansion of the gauge-field part of the supergravity action 
about y4^j = is equivalent to that for a collection of non-interacting (abelian) gauge 



potentials, 



Sa 



(Fx 



M) 



(2.25) 



where a convenient covariant gauge-fixing Lagrangian has been added to remove the 
gauge freedom SA'^.j = Om^""- The gauge and Lorentz indices are raised and lowered 
using the metrics 5ai3 and timn, respectively. 

This then leads to the following momentum-space propagator for A"^ 



Vmn + (6 



kivfk 



(2.26) 



The coupling of any such a 6D gauge boson to a brane again depends on the 
microscopic details. If so, then we expect the linearized coupling to have the general 
form 

5int = Jd"^^ J^'AIj , (2.27) 

where we take the current, J*^, to include the appropriate 6D gauge coupling, g. 

One should also be alive to more complicated possibilities, however, such as if an 
abelian bulk gauge field, F^^^, were to mix with an abelian brane gauge field, V^, 
as in 

S'^, = k^p j d'^ F^,V^^^ , (2.28) 

where k^p represent a matrix of mixing parameters. Phenomenological constraints 
can be quite strong for such interactions, because if they involve the electromagnetic 



field they can produce effects which are easily detected even if very small [|4 



2.5 2-Form Gauge Potentials 

A field which is ubiquitous to higher- dimensional supergravity theories is the skew- 
symmetric gauge potential, Bmn, since this field is usually related to the graviton 
by super symmetry. The quadratic action for this field in the absence of background 
field strengths is 



- / Fx 



- / d^x 



1 

12 
1 
4 



^ r^MNP 



+ 



+ 



1 



{OmB 



MN\2 



2 VCa 



MN\'^ 



1 (c?m5"'^) 



(2.29) 
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where Gmnp = QmSnp + d]^BpM + dpBMN and the second term is the gauge-fixing 
Lagrangian which is added to remove the gauge freedom 6Bmn = dM-^N — dN-^M, for 
arbitrary Am- The corresponding momentum-space propagator for Bmn becomes 



1 



+ 



{VmpVnq - VmqVnp) 
- 1 



(2.30) 



A;2 



(VMpkNkQ - VNpkMkq + rij^qkukp - VMqkNkp) 



The field Bmn can couple to branes, although these couplings often vanish in the 
absence of background field strengths in the bulk or on the brane. An example is the 
Bmn coupling to Z^-branes, which is obtained (in the absence of Ramond-Ramond 
fields) by substituting Gmn — ^ Qmn + i^Bmn + ko'Va/at inside 7^jy in eq. ( |2.19| ), 
where a' is a dimensionful constant and as before Vmn = QmVn — OnVm is the field 
strength of a brane gauge field. 

Other types of antisymmetric tensors can also arise, which are not related to the 
extra-dimensional metric by supersymmetry. Examples of these are the Ramond- 
Ramond fields which arise within Type IIB supergravity models in 10 dimensions. 
(Higher-rank fields can also arise, although we do not pursue the implications of 
these further here.) 

In general, for slowly varying fields we expect a linearized coupling to have the 
general form 



for some model-dependent current, J^^^ 



MN 



(2.31) 



-J 



NM 



2.6 Gravity 

Finally, we follow standard practice and expand the metric about the background 
using 

gMN = Vmn + 2k Hmn ■ (2.32) 
Expanding the Einstein-Hilbert action to quadratic order in Hmn then leads to 



S 



UmH d^'h + U^hpQ + ap/i^^ dMh - dph'''' doh% 



N 



1 

T2 



d^UMN — '^Qn^ 



(2.33) 
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where h = rj^^hMN- With this choice the propagator in momentum space (in 6 
dimensions) is 



k'^ + p 



,2 



^ {VmpVnq + VmqVnp) - ^VmnVpq (2-34) 



(^^2 — 1) 

"^ 2(P + '^^^■^^^^^'3 + VMqkNkp + rj^pkukq + riNqkukp) 

The interaction between these fields and matter is the usual gravitation one, 
with Hmn coupling to the linearized stress tensor: 

jd„ rnMN 



S,^, = -^^ j d^x T^'^'huN • (2.35) 

3. Interactions Amongst Non-relativistic Sources 

In this section we compute the interaction energy which is generated by the linearized 
tree-level exchange of 6D bulk fields between two classical charge distributions. For 
practical applications we focus on the case where the sources sit on the same brane, 
and so are separated only in the 3 visible spatial directions and not separated at all 
in the 2 compact dimensions. 

Consider now the exchange of bulk fields from classical sources described by the 
couphngs of eq. (|]20|), (pl), {^^ , ( p3T| ) and ( g^Sl ). Such an exchange leads to 



the momentum-space interaction potential 

V{k) = ^'[p,{k)D%k)p2{-k) + k''j,Uk)D%{k)k''j2%{-k) 

^MN,PQ 



+Tf^(fc)D,%pQ(A:)r,^«(-A:) , (3.1) 



where k^ = {k^,p} represents the 6-momentum transfer which occurs due to the 
bulk-field exchange. Although all indices written here are 6-dimensional in practice 
the macroscopic sources on the brane have currents which are parallel to the branes, 
which we take also to be the large 4 dimensions x'^. 

For precision tests of gravity our interest is in the elastic scattering of macroscopic 
non-relativistic sources, and so the momentum transfer which is relevant in the large 
four dimensions is k^ = (0,fc). We also keep the extra-dimensional momentum 
transfer, p, arbitrary because here we know the brane states are localized at the 
brane positions, and so the brane-bulk couplings in themselves break momentum 
conservation in these directions. 
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3.1 Point Sources 

In order to apply the above expression to real measurements of forces we must ask 
what the source currents, p, jlj, J^^, Jmn and Tmn, would be for a brane-bound 
source. Although we might expect the answer in general to be model-dependent, a 
considerable amount can be learned from symmetries in the case when the interaction 
is non-relativistic and elastic and the source is labelled simply by its 4-momentum 
and spin. (For instance, in what follows we apply these results to the interactions of 
electrons or protons or any other non-relativistic object whose internal structure is 
not disturbed by the bulk-field exchange.) In this case currents like = {p, s\J^\p, s) 
depend only on the particle 4-momentum and the spin pseudo-vector s^, where 
p'^Pn = —m?, p^s^ = and s^s^ = +1, and so parity and Lorentz invariance deter- 
mines their form up to an overall scalar normalization.^ 

For instance, in such a case the single-particle matrix element of a 4-dimensional 
vector, axial vector, skew-tensor, skew axial tensor and symmetric tensor current 
must have the form 

Jfi = (lP^t, J^ = (lASi,, 

and T^^ = Ap^^p^ + Bs^^s^ , (3.2) 

where g, qa-, c, ca, A and B are unknown constants whose form is not determined 
by 4D Lorentz invariance. 

Furthermore, for static interaction energies our interest is in negligible velocities, 
so we may specialize to the rest frame for which p^ ~ (m, 0) and = (0, s), where 
s is the unit vector pointing in the direction in which the particle's spin angular 
momentum is measured. In this case the m dependence of the various constants 
may be determined on dimensional grounds, leaving only a dimensionless number 
undetermined. Additionally, since all the currents J^, J^,^ and give the same 
kind of contributions to the interaction potential, we only keep J^^, for simplicity. 

Using these arguments we find that the quantities which appear in eq. ( |3.1| ) are 
approximately constant for the range of momentum transfers of interest, X{k) ~ A(0) 

^We impose parity invariance on the matrix elements since the sources of interest are bound 
states of parity-conserving interactions like electromagnetism, the strong force or gravity. 
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etc., and thus 

A(0) = am , 

•^a (0) = gaQaP^/m ~ QaQa K , 

JA'-(O) = I ip'^s'^ - p-^sn = ^ ms^ {6i; 5r - 5o^5f ) , 

T^-(O) = m (^'^o^'^o + ^ -^f ^j) • (3-3) 
Now a, Qa, g^, 6 and c are dimensionless constants, Qa is the 6D gauge couphng 
constant (having dimensions of inverse mass) and we have used the definition of 
mass to fix the normahzation of m in the coefficient A defined in eq. ( p.2|) . 

Inserting these expressions into the interaction potential of eq. ( |3.1| ) (and spe- 
ciahzing for simphcity to a single gauge charge, q) then gives 

\/(F,p) = !^K'm,m2 (x + l)- 9\q2^ , (3.4) 



where 



X = aia2 - CiC2(si -52) + ^ {h + 62) + hh 



(3.5) 

and k'^ = {k^,k). In this expression the first term (0102) arises due to dilaton 
exchange, the last term (giq'2) from the exchange of bulk gauge fields, the spin- 
dependent term (C1C2) from the exchange of the Kalb-Ramond 2-form potential, and 
the (^1,2) terms pluss the penultimate term (3/4) follows from 6D graviton exchange. 
Derivative couplings preclude Goldstone boson exchange from contributing in the 
static limit. 

Taking the limit ^ corresponding to elastic (static) interactions and trans- 
forming to position space we then find 
1 r (fk 



ikr+ip-y 

1 f o 3\ o ) ^ r d^k (^f'-^+ipy 



AurA 



(3.6) 

where r = \r\. Taking two sources which share the same brane we take y = 0, and 
see that the resulting V{r) becomes the expected sum over Yukawa potentials, one 
for each KK mode. 
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For sufficiently large r only the zero mode p = contributes, giving the standard 
1/r result 



^ ^ + 7 + , (3.7) 

^ ' AnrA V 4y ArcrA 

and from this we read off the four- dimensional Newton constant, Gat, by focussing 
on the 6-independent gravitational term (last term in the parenthesis): 

Gn = ^. (3.8) 
16nA 

Similarly, the effective gauge coupling constant, e, is given by 

a = ^ = -l-. (3.9) 

Notice that although eq. ( pl8| ) pre-multiplies just the gravitational contribution 
(so defined because it includes only the contribution of the massless modes of the 
extra-dimensional metric) it differs from what would have been obtained — namely 
Gn = k^/{87tA) — using only the 4D graviton. This difference arises because of the 
contribution in eq. ( p.7|) of other 6D metric polarizations besides the 4D graviton. 

We see here another reflection of the fact that long-distance gravity also deviates 
from General Relativity in these extra-dimensional models, because of the existence 
within them of massless KK metric modes which can mediate long-range interactions. 
Although these same modes also exist in non-supersymmetric models with large extra 
dimensions, it is important to notice that in the non-supersymmetric case quantum 
corrections lift their masses to the generic KK mass scale, thereby removing any 
such long-distance signature. The same need not be true within SLED, since in 
this case the suppression of the quantum contributions to the vacuum energy also 
suppress the corrections to the light-scalar masses, leaving these modes naturally 
light enough to have observable effects over long distances [^]. In this sense long- 
distance modifications of General Relativity potentially distinguish supersymmetric 
from non-supersymmetric large-extra-dimensional models, although the strength of 
the couplings of these modes depends on the details of the extra dimensions and how 
they have evolved over cosmological timescales. 

Although we find in the present instance more than one such massless KK mode, 
this is an artefact of our having compactified from 6D to 4D on a torus. Generic 
compactifications need have only two light scalar fields, corresponding to the dilaton 
arising from the approximate scale invariance of 6D supergravities together with its 
axionic partner under 4D supersymmetry. It is these fields which are ultimately 
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responsible for the time- dependence of the Dark Energy density which SLED models 
predict |T9|, |g. 

Another torus-specific feature of the above results is the fact that the KK spec- 
trum is identical for bulk fields having different spins. For general extra-dimensional 
geometries these spectra can differ, leading to a more complicated interplay of the 
strengths of different fields as a function of source separation. 

With these definitions we finally obtain the following interaction potential energy 
for two point sources situated at the same point in the extra dimensions but separated 
by a distance r in the large 3 space dimensions. 

Vir) = + + ^} E^"'"" • (3-10) 

For general separations we perform the sum in eq. (|3.10|) numerically, with results 
plotted as a function of r in fig. (|l|). As is clear from this figure, for large r the 
interaction potential varies as 1/r, as appropriate for 3 spatial dimensions, while for 
small r the sum over many KK modes combines to give a result which varies as 
as is appropriate for forces in 5 spatial dimensions. The various curves in this figure 
show how the results depend on two of the moduli of the torus: the ratio ri/r2 of 
the lengths of its two sides, as well as the 'twist' angle 6. For very asymmetric tori, 
with one direction much longer than the other {i.e. small ri/r2 or small 6), the figure 
also shows that there can be an intermediate regime for which there are effectively 4 
spatial dimensions, and so for which the force law varies as 

As the figures make clear, the crossover between the 1/r and 1/r^ or 1/r^ power- 
law behaviour occurs over quite a short range of r. This complicates the interpreta- 
tion of experiments, for which deviations to Newton's Law are often parameterized 
using a single Yukawa-type exponential, 

r(r) ~ ^ (1 + ae-'^) , (3.11) 

with the parameters (a and b) found by fitting to the data. Naively, one might expect 
this form to provide a good description, with the parameters found corresponding to 
the properties of the first KK mode. In fact, as shown in fig. (^, the rapid crossover 
to a new power law implies the best fit is usually quite different from what would 
be inferred from the first KK mode, and also depends strongly on the interval of r's 
used in the fit. This emphasizes that some care is required when inferring the form 
of the interaction potential directly from observational data. 



17 




Figure 1: The potential in eq. (3.10) calculated numerically for different choices of the 
moduli 9 and r] = rijr^ (where is defined to be the larger of the two toroidal radii), with 
Vo/(27r) = —(Gt^m\m2lr\) (l + ^) + (pLq\q2lr\). The large- and small-r regimes where 
y(r) ~ 1/r and respectively is clearly seen, as is an intermediate regime for which 
V{r) ~ for very asymmetric torii (ry = 0.01 (left figure) and = 0.01 (right figure)). 
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Figure 2: The exact potential for rj = 1, 9 = 7r/2 (solid line) is fitted to an expression 
with a single exponential, of the form (1 + ae~'"')/r. The best fit (dashed line), obtained 
from a least squares fit to the logarithm of V/Vq, is (1 + 29.2481 e~^'^^'^'^^'')/r, quite far 
from the "correct" expression (1 + 4e~'~)/r from the first KK mode (dotted line). Here r 
is written in units of ri/(27r), and Vq is defined in fig. (|^). 



A final feature of these expressions which is worth mention is the overall sign of 
the deviation from Newton's Law. Non-supersymmetric 6D models involving only 
graviton exchange generically predict that the strength of the interaction potential 
increases relative to Newton's law as r decreases, and this prediction is robust because 
it follows only from the attractive nature of graviton exchange. 
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The same need not be true in supersymmetric models because although scalar 
and graviton exchange generate attractive forces between identical particles, vector 
exchange leads to repulsion. In order to find a weaker force than gravity as r decreases 
it is therefore necessary to have the lowest vector KK mode dominate over the relevant 
range of distances, which could in principle be arranged by ensuring that the lightest 
KK mode correspond to a bulk vector field. The above figures show that this does 
not occur when these are compactified on a torus, however, because in this case all 
bulk fields have identical KK spectra and so the vectors do not dominate the scalars 
and gravitons. However this is an toroidal artefact, and a detection of an initially 
negative residual could carry considerable information about the shape of the extra 
dimensions. 

3.2 Extended Sources 

It is useful to generalize the result of eq. ( p.lOp , from point sources to extended 
sources. This is a necessary complication because once the interaction potential 
between point sources is no longer proportional to 1/r it is not true that the potential 
due to a spherical source distribution is the same as that of a point charge. 

In the present instance the interaction energy of two sources having mass den- 
sities pi and p2 whose centres are displaced by a separation r is calculated by inte- 
grating eq. ( p.lO|) , to give 

-j-p\r+r2-ri\ 

Pi(r'i)p2(^'2)-p 



V{r) = -gy^[ SrJ rfV2Pi(ri)p2(r2)^^^^^^^. (3.12) 
„ ii/i Jvo \r + r2-ri\ 



— auJiUJ2 (3.13) 



Here 

G = a.r ( ^ ^ ^ 

where the uji denote the charge-to-mass ratios for the two source distribution, which 
is assumed to be constant. Similarly {X) denotes the relevant average over the source 
densities. 

The density integrals can be calculated explicitly for simple matter distributions, 
some of which we now perform for illustrative purposes. 

Spherical sources 

If we assume that the two sources are spheres with constant density, pi and p2, and 
radii, Ri and i?2, so = AnpiRf/S. We further assume r > _Ri + _R2 so that the 
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two sources do not physically overlap. With these choices the interaction potential 
energy becomes 

g-p|r+r2-ri[ 



V{r) = -gpiP2y2 I [ d\2 

„ Jr^<R^ Jn<Ro 



ri<Ri Jr2<R2 \f + f2 — 'ri\ 

.^^$:5(|p|i^,)5(|p|i?.)e^-^ (3.14) 



where 1121 



^, , 3 (x cosh X — sinh x) ^ ^'^ A^ 

B{x) = — ^ 1 + — + O(x^) if X < 1 

x'^ 10 

3e^' / 1 



2x2 



(^1 - + O (e-^) if x> 1.(3.15) 



For each KK mode the effect of distributing the matter over a sphere is to 
multiply each factor of mass, Mj, by the function B{\p\Ri). Since lim^-^o -S(a;) = 1 
the contribution of any KK mode whose wavelength, A = l/|p|, is longer than the 
sphere's radius is not modified from that of a point source, and so in particular this is 
true for the KK zero modes. Since these wavelengths are at their largest shorter than 
a micron lengths, this means that for the spheres of practical interest the macroscopic 
effects encoded in B are important for all nonzero KK modes. For these it is the 
opposite limit, \p\R ^ 1, which is relevant, in which case 



gM,M2 

V{r) ^ 



9 



y- 



e 



-lp|(r-iJi-i?2) 



(3.16) 



The sum over p is exponentially small unless the spheres are separated by micron- 
sized distances of order the longest KK wavelength, r — i?i — ^ A ~ max[l/|p|], 
and for separations this large it is 0{\^ j 

The sum in eq. ( ^.14| ) is performed numerically and compared to the point source 
potential in fig. (^), for the simplest possible torus with r\ = r2 and 9 = 7r/2. The 
ratio between the two potentials is always larger than 1 since -B(x) > 1. We can find 
an upper bound to this ratio when the two spheres of equal radius touch (r = 2R) 
by replacing the sums with integrals: 



Epi?^(|p|i?)e-^l^l^ 
This estimate is independent of the moduli of the torus. 



^'.^ </ 4x52(x)e-2^c/x ^ 1.635 532. (3.17) 

l^p ^ Jo 
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Figure 3: The figure to the left shows the ratio between the potential of spherical sources, 
eq. ( 3.14| ), and the corresponding point source potential. The radius is, when going from 



left to right and then down, R = 0.05, 0.1, 0.2, 0.4, 1, 1.5 and 2. The figure to the right 
shows the relative contribution of the KK sum in ( p.l8| ), for the potential between the 
Earth and a small sphere. The size of the torus is here chosen as ri = 10~^ m. For both 
figures all lengths are written in units of ri/(27r), and the moduli of the torus are chosen 
as ri = r2 and 9 = 7r/2. 

Semi-Infinite Bulk Sources 

Of important practical interest is the case where a very small spherical object inter- 
acts with a very large one (say, the Earth). The potential energy in this case may be 
obtained from the above result for spheres by taking the limit where the radius of one 
of the spheres becomes very large. We are led in this way to the interaction energy 
between a sphere of mass M and radius R, whose centre is displaced a distance z 
from the edge of the semi- infinite source. 

Denoting the density and radius of the large sphere by p and L, we define its 
mass-per-unit-area as a = M/A = (47rpL^/3)/(47rL^) = pL/3. Dropping an overall 
additive constant and neglecting sub-dominant powers of z/L we find the interaction 
energy to be 



V{z) = M 



gz-2ngpJ2mP\R)^e-^^\^ 



pi 



(3.18) 



where g = AttQcx. The first term in this result describes the constant (Galilean) 
acceleration towards the large source due to the exchange of massless KK modes, and 
since this acceleration 'sees' the entire source its strength depends on the mass/area 
ratio a. 
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By contrast, the second term represents an attraction of the smaller sphere to- 
wards the nearest piece of the larger sphere due to the exchange of massive KK 
states, which is only sensitive to the amount of matter within a range A of the edge 
of the small sphere due to the finite range of these interactions. Because of this finite 
range this part of the interaction is controlled by p rather than o", and is exponen- 
tially small unless the sphere's edge is within A {i.e. closer than a micron or so) from 
the edge of the semi-infinite source (see fig. (^). Because this attraction is only to 
nearby source material, it is also sensitive to any local deviations of the large source 
from spherical symmetry. 

The apparent divergence in eq. (|3.18| ) as 2; ^ (due to the divergence of the 
sum) is illusory so long as R is nonzero. This is because the geometry requires z > R, 
with the sphere touching the slab only when z = R. The result is finite in this limit, 
as may also be seen explicitly from its parent formula, eq. ( 3.16 ). Eq. ( 3.18 ) would 
diverge in the limit i? — 0, however, since in this case B = 1 and z can approach 
zero. This represents a divergence in the potential of a single spherical which only 
arises at the source's edge, and as such it may be absorbed into the coefficient of an 
effective interaction which is localized on the boundary of the source. 



In the limit of practical interest, \p\R ^ 1, we may use in eq. (|3.18|) the asymp- 
totic form for B{x) to obtain the simple formula 



V{z) = M 



gz 



3nQp 



,-\p\d 



(3.19) 



where d 



■ R denotes the separation between the nearest edges of the two sources. 



4. Conclusions 

The calculations in this paper allow the following conclusions to be drawn concerning 
the implications of SLED models for tests of gravity at short distances. 

• SLED models predict deviations from the inverse- square law at micron dis- 
tance scales, giving instead a position-dependence which describes a crossover 
from 1/r^ to 1/r^ over a fairly short range of distances. These predictions 
are not well- described by perturbing Newton's Law with a single Yukawa-type 
exponential potential. This situation is similar to what is obtained in non- 
super symmetric models with sub-millimeter extra dimensions, with the impor- 
tant difference that within SLED the scale of the crossover is related to the 
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observed Dark Energy density and so cannot be adjusted to smaller values if no 
deviations from Newton's Law are observed. It also differs from the Yukawa- 
type deviations which are likely to be produced by recent string models having 



a string scale in the TeV region with six large dimensions ||26 



• The precise shape of the deviations from Newton's Law in the crossover region 
are likely to depend on the details of the Kaluza-Klein spectrum, and so also 
on the precise shape of the extra-dimensional geometry. This is already visible 
in the calculations presented here, through the dependence of the forces on 
the various shape moduli of the extra-dimensional torus. Calculations are 
underway to determine more precisely how this shape depends on these details 
for more general geometries than the tori considered here. 

• SLED models also differ from their non-supersymmetric counterparts in the 
number and kinds of fields which supersymmetry dictates must populate the 
extra dimensions. In particular, since bulk vector exchange mediates a repulsive 
force, it is in principle possible to obtain an interaction energy which is weaker 
than Newton's law at short distances if the lightest KK mode should be for 
the vector fields rather than the tensors or scalars. So far as we know this 
same possibility cannot arise within the purely-gravitational extra-dimensional 
models considered heretofore. This type of short-distance repulsive interaction 
turns out never to dominate for the simple toroidal compactifications considered 
here, but could be possible for more complicated extra-dimensional geometries. 

• The presence of several bulk fields in SLED models can also generically intro- 
duce a composition dependence to the observed inter-particle force. Unfortu- 
nately, the detailed nature of this composition-dependence is difficult to predict 
without knowing more of the details of our extra-dimensional surroundings. 

All of these conclusions underline the fact that the SLED proposal is unique in 
that it provides a very precise framework within which a fundamental connection 
is made between the observed Dark Energy density and other kinds of observable 
physics. In it deviations of gravity at sub-millimeter distances are also tied to ob- 
servable signals at the LHC and possibly to scalar-tensor deviations from General 
Relativity over both astrophysical and cosmological distance scales. We hope that 
experimenters take up the challenge of searching for these phenomena. 
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